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Abstract. I present here a new Glauber-inspired approach to derive normal cold 
nuclear matter effects on J/ip production in Au-Au collisions. In an as much as possible 
Cn| ' model-independent way, it extrapolates the centrality dependent yields from d-Au to 

Au-Au collisions. I then compare to the new Au-Au measurements shown by the 
| PHENIX experiment in this conference. In the most central collisions, J/i/> survival 

£S| ' probabilities beyond cold nuclear matter effects are found to be 44 ± 23% at forward 

Q : ( 9 = 1.7) md 25 ± 12% at Md(v = 0) rapidities. 
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1. Cold nuclear matter effects 



J I ip production is sensitive to nuclear matter and subtracting the normal suppression is 
a key point in any interpretation of quarkonia production in heavy ion collisions. Data 
from p-A or d-A collisions provide good insights into cold nuclear matter effects. 

At SPS energies, a single J / '-^-nucleon absorption cross-section of o 3 J^ = 4.18 ± 
0.35 mb can describe the observed J/ip suppression in a wide range of p-A systems, as 
well as in S-U and peripheral Pb-Pb or In- In collisions [U [2]. However, according to 
various theoretical predictions (such as [3]), shadowing (or rather anti-shadowing) could 
play a role at the SPS regime, namely for momentum fraction x of the order of 10 _1 . The 
SPS experiments can hardly address this issue because of their limited rapidity (or x) 
range. However, a large rapidity asymmetry observed in J/ip yields in p-A collisions 
provides a hint that a nuclear effect is playing a role beyond simple absorption [I] . 

At RHIC energies, the PHENIX experiment has measured [5] J/tp nuclear 
modification factors in d-Au collisions (RdAu)- It covers a wide rapidity range, from 
—2.2 to 2.4, corresponding to momentum fractions x of gluons from 10~ 3 to 10" 1 . It 
also includes the first ever centrality dependence measurement of J/ip production in 
p-A like collisions. According to [Bj, both the rapidity and centrality dependencies can 
be reproduced in terms of inhomogeneous shadowing added to a moderate (not larger 
than 3 mb) normal nuclear absorption. In [7j, the same author derives Au-Au nuclear 
modification factors (RauAu) , assuming various shadowing schemes and absorption cross- 
sections. Another attempt to derive a reference from d-Au can be found in [8] where 
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Figure 1. Phenomenological fits to RdAu versus impact parameter. Left) linear 
fits and one standard deviation variations. Right) other fits for systematic studies: 
Aexp(6/B), I- (I- A) exp(-6/-B), A + (1 - A) tanh(6/B) or exp(-AVi? 2 - b 2 ). 

the authors compute dissociation cross-sections from the centrality dependence of RdAu- 
They assume that all nuclear effects are proportional to exp(— p (j diss L), p being the 
normal nuclear density and L the average length of nuclear matter seen by the J/i/j. 
However, nuclear effects such as shadowing do not follow this centrality dependence. 

2. From d-Au to Au-Au at RHIC 

I propose here an alternate method, with the concerns to be as much data-driven as 
possible and to derive realistic uncertainties. To do so, I first perform phenomenological 
fits of the modification factors RdAu{y, b) as a function of the impact parameter b (given 
by a Glauber model), for the three rapidities y probed by PHENIX [5J. The simplest 
working fits are linear (figure [D left). A cut-off at RdAu = 1 is applied at y = and 
y = +1.7 for high impact parameters to reflect the fact that grazing Au-Au collisions 
are expected to behave like (binary-scaled) p-p collisions. This will prevent RauAu from 
reaching unphysically high values (well above unity) at high impact parameters. 

I then run a Au-Au Glauber model [9] simulation. For each Au-Au collision 
occurring at a given impact parameter b Au Au, the N coU elementary nucleon-nucleon 
collisions are randomly distributed (following Woods-Saxon nuclear densities) providing 
the locations b\ and b\ of each collision i, relative to the center of nucleus 1 and nucleus 2. 
For the considered Au-Au collision, the predicted nuclear modification factor RauAu is 
given by the following summation over the elementary collisions: 

N co u 

RAuAu(y,b AuAu ) = ( R dAu(-y,b\) x R d Au(+y,b l 2 ))/N co u. (1) 

1=1 

This formula assumes that a J/ip produced in a Au-Au collision at a given rapidity 
y suffers the product of the nuclear effects that were observed in d-Au at this rapidity, 
by the ones of the opposite rapidity (equivalent to a Au-d collision). This assumption is 
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Figure 2. RauAu versus the number of participants N part at y — 1.7 (left) and y = 
(right). Theoretical curves are nuclear effects predictions from Vogt [7], solid and 
dashed lines being for a a i, s ~ 1 and 3 mb. The circles within the shaded bands show 
the prediction from the model presented here. Data points are from PHENIX [ID] . 

at least correct for the only two effects considered so far to explain RHIC data, namely 
shadowing and nuclear absorption, since quarkonia production is then proportional to 
the parton distribution functions (pdf) in each nucleus, while the average length is the 
sum of the length in each nucleus: RauAu oc pdfi x pdf 2 x exp(— pa{L\ + L 2 )). 

The method also allows to fully propagate the errors from the d-Au measurements, 
both statistical and uncorrelated systematic, to the Au-Au prediction, by varying the 
RdAuip) uncorrelated fit parameters by one standard deviation, as shown by dashed lines 
on figure [U left. In addition to that, two sources of systematic errors are considered. 
First, the fit function is not theoretically known and I tried various shapes, as displayed 
on figure [1] right. In particular, this helps to cope up with the lack of sampling 
in the low impact parameter region. The obtained RauAu vary by up to 5 or 12% 
at forward or midrapidity, asymmetrically and depending on centrality. Second, the 
Glauber parameters (total cross section, Woods-Saxon parameters,...) are varied within 
allowed limits,which results in RauAu modifications of 2 or 4% at forward or midrapidity. 

Since the approach is based on a Glauber calculation, it is easy to average the 
predicted modification factors on experimental centrality classes. This is done on figure [2] 
where predictions are given for the Au-Au PHENIX centrality classes from [10] and 
compared to the measured RauAu- At both rapidities, they are comparable to Vogt's 
predictions and incompatible with the measured RauAu, pointing out that J ftp are 
suppressed beyond cold nuclear matter effects, namely anomaly suppressed. 

This method allows proper propagations of the experimental uncertainties and 
selections of centrality classes. Its main disadvantage lays in the fact that, to predict 
a given Raa, it requires p-A (or d-A) collisions at the same energy, with the same A 
nucleus and for both positive and negative rapidities. In particular, these are not yet 
available for Cu-Cu collisions at RHIC, and not foreseen at LHC. 
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Figure 3. J/ip survival probabilities beyond cold nuclear matter ar forward (squares, 
y = 1.7) and mid (circles, y — 0) rapidities. 



3. Survival probabilities beyond cold nuclear matter 

On figure [3j 1 divide the observed suppression by the predicted nuclear effects. The boxes 
stand for the Au-Au correlated errors, enhanced by the (dominant) cold nuclear matter 
uncertainties, while the bar is for the Au-Au statistical and uncorrelated systematic 
uncertainties. Additional global uncertainties of 30 or 35% for forward or midrapidity 
need to be considered, mostly due to the fact that the p-p references used in RdAu and 
RauAu are different and do not cancel in the ratio. A reanalysis based on renormalized 
RdAu to the p-p reference used in RauAu will substantially reduce this to ~10%. 

Taking all the uncertainties into account, the most central survival probabilities are 
as low as 44 ± 23% at forward (y = 1.7) and 25 ± 12% at mid (y = 0) rapidities. 

It's my pleasure to thank Jean Gosset and Klaus Reygers for valuable help to get 
this work started. 
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